Prom electronic structure to catalytic activity: A single descriptor for adsorption and 

reactivity on transition-metal carbides 
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Adsorption and catalytic properties of the polar (111) surface of transition-metal carbides (TMC's) 
are investigated by density-functional theory. Atomic and molecular adsorption are rationalized with 
the concerted-coupling model, in which two types of TMC surface resonances (SR's) play key roles. 
The transition-metal derived SR is found to be a single measurable descriptor for the adsorption 
processes, implying that the Br0nsted-Evans-Polanyi relation and scaling relations apply. This gives 
a picture with implications for ligand and vacancy effects and which has a potential for a broad 
screening procedure for heterogeneous catalysts. 
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Finding new simple, efficient, and cheap catalysts is 
one way to solve some of today's global environmental 
challenges 0], as illustrated by a recent way to harvest 
solar energy A goal in catalysis research is to de- 
sign and tune the activity and selectivity of catalysts by 
controlling their structural properties at the atomic level. 
This calls for an identification of key concepts, which can 
be done from fundamental theory, in particular density- 
functional theory (DFT), as shown for transition- metal 
(TM) surfaces [3]. 

For such surfaces, the important role of the TM d elec- 
trons for chemisorption and catalysis was early antici- 
pated yi, estimated [5|], and articulated in the d-band 
model [sl, 0] , which establishes relations between atomic 
structure and activity for TM catalysts The d-band 
model explains and predicts a variety of TM properties, 
for example, adsorption energies on TM alloys, strengths 
of bonds at steps and terraces, and transition-state en- 
ergies [1, 01, by correlating them to the energy of the 
d-band center. The success of the d-band model has led 
to the introduction of simple descriptors that are able 
to rationalize experimental data Examples of such 
descriptors, e.g., for the water gas shift reaction, are the 
adsorption energies of oxygen and carbon monoxide on 
TM surfaces. Today, design of new TM catalysts by com- 
putational screening is a realistic approach Q . 

A need for catalysts, a scarcity of precious catalyst ma- 
terials, and a general curiosity spur the interest in cat- 
alyst materials beyond the TM's. Statements like "for 
several types of reactions, such as hydrogenation reac- 
tions, catalytic activities of carbides and nitrides were 
approaching or surpassing those of noble metals" can be 
found in the literature [lO|. Thus, a theoretical approach 
extended to more complex materials than TM's is of in- 
terest. In a study of H adsorption on TM-terminated 
(111) surfaces of TM carbides (TMC'sV a deviation from 
the d-band model has been observed [ll| . A recent study 



of TM oxides, nitrides, and sulfides stresses similarities in 
scaling behavior between molecular and atomic adsorp- 
tion energies but calls for "a suitably modified d-band 
model" [13 ■ There exist suggestions of a bulk-derived 
descriptor for the reactivity of the TMC's ,13;], however, 
our work goes beyond and finds a surface-derived descrip- 
tor based on electronic-structure calculations. 

Electronic surface states and resonances (SR's), known 
to appear at steps and defects, play an important role for 
the catalytic activity of a system. The ideal and stable 
TMC(IOO) faces are common objects of study in the liter- 
ature [ij] and do not show any presence of SR's. Recently 
we have shown that the understanding of atomic adsorp- 
tion on TiC and TiN calls for substantial steps beyond 
the original d-band model [H, H, [i3| • On TiC(lll) and 
TiN(lll), the key actors are identified to be two types 
of surface resonances (SR's), derived from Ti (TiSR) 
and C(N) states (CSR/NSR), respectively. The adsorp- 
tion mechanism is explained within a concerted-coupling 
model (CCM), in which there are two types of interac- 
tions involved, one between the adsorbate and the TiSR 
and one between the adsorbate and the CSR's (NSR's) 
Til IgI . 17 1 . The presence of SR's motivates a focus on 



the metastable and polar TMC (111) surfaces, which may, 
for example, serve as prototypes for steps and defects. 

This Letter demonstrates that an approach that in- 
volves DFT calculations, trend studies, analyses of com- 
puted results in electron-structural terms, formulation of 
a model, and utilization of model predictions 15, 1^, 17 1 
is useful to reach the above stated goal. In terms of the 
CCM, we identify a single descriptor, the mean energy 
of the TM-localized SR (TMSR), that accounts for three 
important processes on the considered TMC surfaces: ad- 
sorption, dissociation, and catalytic activity. The results 
should be general but are here illustrated for one surface 
reaction, the ammonia synthesis. As a consequence of the 
single descriptor, calculable and experimentally available 
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for some substrates [l8j, scaling relations and Br0nsted- 
Evans-Polanyi (BEP) relations are shown to apply, in 
the same way as for TM surfaces 0. Finally, the gener- 
ality and practical usefulness of the conceptual picture is 
demonstrated by applications to ligand and defect effects, 
which also indicate possibilities for further refinements. 
Taken together, these results provide a framework for a 
systematic analysis of the catalytic activity of the TMC's. 

Trend studies are conducted for the catalysis-relevant 
adsorbates H,N, O, and NH^: {x = 1, 2, 3), with respect to 
the TM component of the substrate TMC, spanning three 
different periods and four different groups in the periodic 
table (TM = Sc, Ti, V, Zr, Nb, Mo, Ta, and W) O. The 
first-principles study is performed within the density- 
functional theory (DFT), as implemented in DACAPO 



[2Ct |2l|. The computational procedures are favorably 
tested and compared with literature data [2^ . From the 
DFT calculations, the adsorption is described in terms 
of adsorption energies (-Bads) and atom-projected local 
densities of states (LDOS). 

Calculated total DOS and LDOS's for the TM- 
terminated TMC(lll) surfaces (illustrated for some rep- 
resentative cases in Fig. [T]) reveal the existence of SR's 
on all surfaces. There is a TMSR in the vicinity of the 
Fermi energy (Ep) and there are several carbon-localized 
SR's (GSR's) deeper down in the valence band [2J]. Note 
that ScC(lll) has no filled metallic d SR states. 

The TMSR's are characterized by studying the differ- 
ence in bulk and surface DOS's [24|, where they appear 
as positive peaks owing to the build up of states at the 
surface. A parameter ecCM is defined as the mean energy 
(center of gravity) of the TMSR, in close analogy with the 
d-band center in the d-band model. The integration 
is performed over the energy range of the positive TMSR 
peak. The value of eccM decreases as the group num- 
ber of the TM component increases (Fig. [H), as expected 
from the filling of the TM d states 0. 

Analysis of the DOS before and after adsorption (illus- 
trated for 0/VC in the lower panel of Fig. [1]) shows that 
the TMSR is depleted upon adsorption. This supports 
the CCM assumption that only the surface-localized part 
of the d-band spectrum is relevant for the adsorption 
and that the TMSR's are key actors in the adsorbate- 
surface interaction [24|. Hence, eccM is a descriptor for 
the atomic adsorption on TMC's. An exception is ScC, 
with its empty TMSR and consequently an adsorption 
mainly arising from the interaction with the CSR's 24 1 . 
This exception confirms the presence of two types of in 
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teractions, as formulated in the CCM 

The role of £ccm as a descriptor is confirmed by ex- 
tensive DFT calculations, which yield a linear correlation 
between the i?ads and £ccm values for each of the studied 
adsorbates (Fig. [2|). The iJads value decreases {i.e., the 
adsorption gets stronger) as £ccm increases, with a gra- 
dient that varies strongly between different adsorbates. 
The slopes of the i?ads vs. eccM correlations are steeper 




FIG. 1: Upper panel: Calculated atom-projected and total 
(111) surface DOS's for a group of TMC's. The dashed di- 
agonal line connects the mean energies eccM of the TMSR's 
and illustrates the downward shift of the TMSR energy along 
the TMC series as the group number of the TM component 
increases. Lower panel: Difference in surface DOS induced by 
O adsorption on VC(lll). 



for the TMC(lll) surfaces than for TM's From a de- 
sign perspective, this could be a useful property, allowing 
larger effects from small changes. Again, ScC does not 
follow the trend of the other TMC's, reflecting the dif- 
ferent nature of its interaction with the adsorbate. 

To further test that proper key concepts (TMSR and 
£ccm) have been identified, ligand effects are studied. 
Manipulating the local environment of the surface TM 
atom changes ecCM, thereby changing adsorption, activa- 
tion, and other energies, which opens up for optimization 
of the catalytically active site. In our illustrative exam- 
ple, O adsorption on a TiC(lll) surface with one, two, 
or three surface Ti atoms next to the fee adsorption site 
replaced by V atoms, the calculated adsorption strength 
is successively reduced at the same time as the local value 
of ecCM is shifted down. This is expected in the CCM 
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FIG. 2: Adsorption energies of H, N, O, and NH^ {x = 1, 2, 3) 
on the considered TMC(lll) surfaces vs. the descriptor eccM 
(the TMSR center of mass). The adsorption energy is defined 
as E^^g = Etmc+a — Etmc — Ea for each adsorbate A. (ScC is 
not included in the linear regressions) . Also included (dashed 
line) are the transition-state (TS) energies for N2 dissociation, 
calculated with a fixed bond-length method. 



(exchange of Ti with V shifts d levels down in energy) 
and monitored by analyzing our calculated LDOS's also 
in these cases. 

However, there is a deviation from the linear -Bads vs. 
£CCM relation [25|. For example, the i?ads value for the 
system with three neighboring V atoms on the TiC (111) 
surface is similar to that for pure VC(lll), although their 
£ccM values differ. This ligand calculation points at fur- 
ther possible refinements of the descriptor. An orbital 
projected DOS shows that the main interaction takes 
place between the adsorbate and the dxz+yz components 
of the TMSR (the z direction being perpendicular to the 
surface). Hence a refined descriptor, taken as the center 
of mass of these levels, could be introduced. However, 
such a detailed descriptor level might not be practical 
from an experimental point of view. 

As ecCM is a descriptor for both atomic and molecular 
adsorption (Fig. [2]), there is a scaling relation between 
the molecular and atomic Bads values (Fig. [3]). Similar 
scaling relations have been found for several molecular 
species on TM's [23| and for OH^ on TM oxides, NH^ 
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FIG. 3: Adsorption energies on TMC(lll) surfaces (ScC is 
not included) for NH^ (x= 1,2,3) vs. N atomic adsorption 
energies (scaling relation). The adsorption energy is defined 
^ -^iTcS^ = -Btmc-i-nh^ — -Etmc — -Enh^. Inset: the linear 
variation of the activation energy for N2 dissociation (-Ets) 
vs. N atomic adsorption energy on TiC, VC, and (5-MoC (BEP 
relation) . 



on TM nitrides, and SH^ on TM sulfides [ij]. These 



scaling relations can be described by a linear equation 
7(^)-^ads + ^' where 7(x) depends only on the 



E. 



AH^ 
ads 



number x of H atoms in the molecule [3|. Our results 
(Fig. [3]) yield values (7nh = 0.62, 7NH2 = 0.33, 7NH3 = 
—0.13) that are close to the predictions 2/3, 1/3, and 
for a; = 1, 2, and 3, respectively. The small deviations 
are likely to be caused by the C species in the substrate, 
which influences the adsorption mechanism directly via 
the GSR's and indirectly via the interaction with the TM 
atom in the compound Jill, [l5, 3 17|. 

In heterogeneous catalysis, linear correlations (BEP re- 
lations) between activation and adsorption energies play 
an important role . Activation energies for N2 dissoci- 
ation on some TMC(lll) surfaces [i^l, calculated as the 
transition-state energy barriers (Bts) by use of a fixed 
bond-length method, are presented in the inset of Fig. [3] 
as a function of the N adsorption energy. They show that 
the BEP relation holds for the TMC's, as expected from 
the linear correlation between the activation energies and 
£CCM (dashed line in Fig. [2|). The result is expected to 
apply generally, as variations in activation and adsorp- 
tion energies are governed by the same basic mechanism, 
which for TM's is the d-band model Si and for TMC's 
the COM [11, 
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The linear dependence of both adsorption and activa- 
tion energies on eccM opens up the possibility to design 
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the TMC substrates using the CCM to suit different cat- 
alytic reactions. This can be illustrated by the ammonia 
synthesis (N2-I-3H2 ^ 2NH3), which is a commonly used 
reaction to develop new concepts and ideas in catalysis 
28 1. By using a slightly modified micro- kinetic model 
29, [s^l the calculated catalytic activities of the TMC's 
are found to order as a volcano curve with respect to the 
single descriptor £ccm 



m 



The descriptor eccM is derived here for TMC's but is 
expected to be directly applicable for understanding pro- 
cesses on other TM compound surfaces (and probably 
even on other surfaces) with SR's similar to the ones on 
the TMC(lll) surface. Such SR's can appear on stable 
surfaces near defects, such as vacancies, steps, and ad- 
sorbed clusters. For example, experimental studies show 
that oxygen atoms replace carbon atoms on TiC(lOO) and 
ZrC(lOO) surfaces when exposed to O2 [3, Our cal- 
culations on an O atom adsorbed in a C vacancy site on 
the TiC(lOO) surface show the presence of a pronounced 
TMSR localized around the C vacancy. The calculated 
adsorption energy (—8.71 eV) is close to the value on the 
TiC(lll) surface (—8.76 eV), to be compared with the 
value (—4.96 eV) on the vacancy-free TiC(lOO) surface, 
where no TMSR is present. This illustrates the general 
importance of TMSR's for surface reactions. 

Taken together, our results provide a framework for a 
systematic analysis of the catalytic activity of the TMC's. 
We have shown that by theoretical means it is possible 
to find a single descriptor, corroborated by experimental 
results on the electronic structure of the TMC surfaces 
[ist . for atomic and molecular adsorption, as well as for 
a simple reaction on the TMC's. The found scaling and 
BEP relations provide simple and efficient connections 
between the (local) surface electronic structure and its 
reactivity, of importance for studying defect {e.g., va- 
cancy) and ligand effects. There are thus implications 
for surfaces, nanosystems, and catalysis, and the intro- 
duced concepts should ramify to other systems and stim- 
ulate the development of similar models for other classes 
of materials. 

Acknowledgments. The calculations were performed at 
HPC2N and NSC via the Swedish National Infrastruc- 
ture for Computing. Support from the Swedish Research 
Council is acknowledged. B. I. Lundqvist gratefully ac- 
knowledges support from the Lundbeck foundation (Den- 
mark) via the Center for Atomic-scale Materials Design. 



* Corresponding author: | alevoj@chalmers.se] 

[1] D. Nocera, Photosynth. Res. 91, 133 (2007). 

[2] D. Nocera, Science 321, 1072 (2008). 

[3] T. Bligaard et ai, in Chemical Bonding at Surfaces and 
Interfaces, edited by A. Nilsson, L. G. M. Pettersson, and 
J. K. N0rskov, (Elsevier, Amsterdam, 2008). 

[4] D. M. Newns, Phys. Rev. 178, 1123 (1969). 



[5] B. I. Lundqvist et al. Surf. Sci. 89, 196 (1979). 
[6] B. Hammer and J. K. N0rskov, Nature 376, 238 (1995). 
[7] B. Hammer and J. K. N0rskov, Adv. Catal. 45, 71 (2000). 
[8] J. Sehested et al. Top. Catal. 45, 9 (2007). 
[9] F. Studt et al. Science 320, 1320 (2008). 
[10] J. G. Chen, Chem. Rev. 96, 1477 (1996); H. H. Hwu and 

J. G. Chen, Chem. Rev. 105, 185 (2005) and references 

therein. 

[11] J. R. Kitchin et al, Catal. Today 105, 66 (2005). 
[12] E. M. Fernandez et al, Angew. Chem. Int. Ed. 47, 4683 
(2008). 

[13] H. Toulhoat and P. Raybaud, J. Catal. 63, 216 (2003). 

[14] F. Vines et al, J. Chem. Phys. 122, 174709 (2005); F. 
Viiies et al, J. Phys. Chem. Ill, 1307 (2007). 

[15] C. Ruberto et al, Phys. Rev. B 75, 235438 (2007). 

[16] A. Vojvodic et al. Surf. Sci. 600, 3619 (2006). 

[17] C. Ruberto et al. Surf. Sci. 600, 1612 (2006); Solid State 
Commun. 141, 48 (2007). 

[18] K. Edamoto et al, Phys. Rev. B 46, 7127 (1992); K. 
Edamoto, et al, Phys. Rev. B 43, 3871 (1991); T. 
Anazawa et al. Surf. Sci. 328, 263 (1995) 

[19] To reduce the complexity and allow for trend studies, all 
considered TMC's are described in NaCl structure, which 
all except WC adopt in either stable or metastable phase. 

[20] S. R. Bahn and K. W. Jacobsen, Comput. Sci. Eng. 4, 
56 (2002); http://wiki.fysik.dtu.dk/ dacapo 

[21] The spin-polarized approximation of the exchange- 
correlation functional according to Perdew-Wang 
(PW91) [2^ is employed. Ultra-soft pseudopotentials 
[33 ] and 400 eV energy cutoff for the one-electron 
orbitals are used. The fe-point sampling of the 3x3x8 
TMC(lll) slabs is performed in a (4,4,1) Monkhorst- 
Pack grid. Repeated slabs are separated by at least 
10.8 A thick vacuum region. During structural relax- 
ations, the two bottom atomic layers are constrained to 
the corresponding bulk geometry. Adsorption is allowed 
on one (the TM terminated) of the two slab surfaces in 
the fee site (the stable site for adsorption on TiC 15]) 
and the electrostatic potential is adjusted accordingly 
[3^ . The influence of higher coverage has been consid- 
ered in Ref. [l^, illustrating the dominance of direct 
adsorbate-substrate forces over lateral interactions. 

[22] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[23] The calculated TMC lattice parameters in the bulk 
(NaCI) phase are a = 4.684, 4.332, 4.164, 4.702, 4.492, 
4.450, 4,479, 4.382 A for ScC, TiC, VC, ZrC, NbC, MoC, 
TaC, WC, respectively. 

[24] A. Vojvodic and C. Ruberto (unpublished). 

[25] For the ligand systems with one, two and three V atoms 
we obtain £ads = -8.09, -7.78, -7.28 eV and eccM = 
0.01,-0.08,-0.11 eV, respectively. 

[26] F. Abild-Pedersen et al, Phys. Rev. Lett. 99, 016105 
(2007). 

[27] Owing to a strongly adsorbed molecular precursor, rel- 
evant activation energies are calculated with respect to 
the molecular adsorption energy. 

[28] A. Hellman et al, J. Phys. Chem. B 110, 17719 (2006); 
M. Boudart, Top. Catal. 1, 405 (1994). 

[29] It is assumed that N2 dissociation is the rate-determining 
step and that all other elementary steps are thus in quasi- 
equilibrium 3013. The turnover frequency (TOP) is given 
by TOF= 2fcPjv2&»&jv2 (1 — 7), where 7 is the approach 
to equilibrium. The BEP and scaling relations are used 



5 



when calculating the TOF as a function of eccM- 
[30] K. Honkala et ai, Science 307, 555 (2005). 
[31] The ammonia synthesis is used just to prove the benefits 
of a single descriptor, not to suggest that TMC's are par- 
ticularly good catalysts for it. Rather, reactions requiring 
a strong molecular precursor, for instance, oxidation of 



methane, are a more interesting group of reactions in a 

real application 
[32] Y. Shirotori et a/., Surf. Sci. 584, 237 (2005). 
[33] D. Vanderbilt, Phys. Rev. B 41, R7892 (1990). 
[34] L. Bengtsson, Phys. Rev. B 59, 12301 (1999). 



